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Abstract 



'■^ , Shot noise and higher current moments Tm ~ (Tr(i^i)^^) are stud- 

PlJ ' ied within the Anderson model of disordered conductors in dimensions 

O ' d = 2, 3 and 4. Here t denotes a transmission matrix. We calculate the 

, >■' , , conductance G = Ti, shot noise S = Ti ~ T2, Fano factor F = S/G and 

ratios Gm = Tm/G, and we show indications that all of them are as good 
order parameters as the conductance itself. In the limit of infinite system 
^ . size, two limiting values of F and Cm are found; the stable one in the 

04 ' metallic regime, and the unstable one, characterizing the critical point for 

d > 2. We present analytical expressions for both limiting values, together 
y^ • with a compact formula for current cumulants at 3D criticality. Our data 

confirm also Nazarov's microscopic theory [PRB 52, 4720 (1995)], as we 
show numerically for special linear combinations of Tm- 



1 Introduction 



d . 

r^ , Much effort has been devoted to the Anderson localization [1]. The transport in 

'"^ ' disordered systems exhibits many universal features, hke universal conductance 

^ I fluctuations (UCF) in the metallic regime, or In G - type distribution in the 

O ■ insulator. 

Theoretical approach concentrated mainly, but not exclusively, on one-dimen- 
sional (ID) models. The type of scaling of the conductance distribution was 
investigated in ID [2] and numerically in 3D [3]. Muttalib et al. [4] analyzed 

}J] ■ and generalized the DMPK equation also beyond quasi-lD and calculated the 

5^ I distribution function of G in all regimes. They tried to apply the generalized 

DMPK also to localized states [5]; in the metallic case it was applied success- 
fully before [6]. Beenakker and Biittiker [7] also studied the metallic transport 
in dimension d > 1 and we will utilize their results for all current moments. 



Besides the UCF, i. e. the second cumulant of G (or var G), numerous 
other exact quantitities where proposed in the metalhc regime (see this Intro 
and Section 2), but only few in the other regimes. The main purpose of this 
work is to calculate a complete set of quantities at 3D metal-insulator transition 
(MIT) exactly. It contributes to theoretical description of the criticality from 
different side than microscopic theory, RMT, non-linear a-model, DMPK, etc. 

Besides the Landauer's conductivity G = Tr T (in units 2e'^/h), and shot 
noise power S — Tr[T — T^], where T = tH, also the universal properties of 
higher order current moments in terms of T^ were observed [8]. The third 
current cumulant was given analytically also in the limit of frequency and tem- 
perature both zero: 

((.=» ^ il') - mm + 2(.)= . (/) (T.-[r(l--^TKl-2T)|) ^ (^ ,^, 

(...) means ensemble averaging. This formula holds for large systems in diffusive 
(metallic) regime and we will reproduce it later. 

Let us denote the higher current moments Tm — (Tr [T*^]) . These quanti- 
ties were analyzed in the works [9] and [10]. The exact solution for normalized 
moments on ID Anderson chain 



Tm r(M - 



1^2 
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Cm - -7^ - T./.,^^9-^./1^9 (2) 



Ti r(M)2r(i)2 

was confirmed by numerical simulations. This formula applies exclusively to 
pure ID. Analogical exact expressions for higher dimensions were given in [7] 
for the metallic transport, Eqs. (22) or (8), and we will add those for 3D critical 
regime. 

Let us make a note on relationship of cumulants and moments. Compact 
expressions for charge (i. e. also current) cumulants <C Q'' S> in diffusive 
regime were derived in [11]. These cumulants can be calculated with help of a 
characteristic function 

OO f^ 

where x is Fourier transformed probability distribution of charges. Assuming 
mutual independence of channels they found 

« Q' » ex J2iTil T)^y-'T\^^^^ (4) 

j 

where the sum over eigenvalues means trace. We get shot noise as the second 
cumulant, proper part of Eq. (1) as the third one, etc. Thus the current 
cumulants and moments are connected in a different way than e. g. those ones 
of conductivity. 



The moments Tm themselves becom.e usually either zero (insulator) or in- 
finite (metal) in the large system limit, except criticality. But we will show, 
that both their ratios and appropriately chosen linear combinations will have 
non-trivial finite values. The limiting Cm ratios have the pleasant property that 
they are independent of boundary conditions. 

The shot noise S = Ti — T2 and Fano factor F = 1 — C2 of Anderson model 
have been studied in [12] for dimensions d = QID, 2 and in [13] for d = 3. QID 
means quasi-lD systems, i. e. stripes of fixed, but large enough width and 
even much larger length. We continue in 2D and 3D with larger systems and 
statistics, and we add both d = 4 and higher current moments M > 3. Under 
simple assumptions on the distribution function of Lyapunov exponents, one 
gets the limiting values for Cm both in diffusive and critical (3D) regimes. We 
will also confirm the validity of Nazarov's theory [14] in diffusive case for higher 
moments and in d = 4; for M =- 2 and d = QID, 2, 3 it was already shown 
in [12] and [13], though less precisely because of smaller system sizes and thus 
larger disturbing ballistic effects. 

There is a great amount of literature on shot noise, e. g. [7], [8]. The Pois- 
sonian value F — > 1 corresponds to the strongly localized (insulating) regime, 
the well known 1/3 suppression F = 1/3 is typical for the metallic regime. If 
d > 2 we get yet another value Fc, which marks the MIT. All higher moments 
will show qualitatively the same behaviour. 

In recent paper [15] even the fourth cumulant, or the combination C3 — C4 
was introduced. Such quantities can be calculated in general. 

The aim of this work is to show that Cm and special linear combinations of 
Tm converge to universal values in the large system size limits; to observe two 
non-trivial limiting Cm values for d > 2; to give them analytically (except 4D 
criticality) and to prove them numerically. 

The paper is organized as follows. In this Intro and in Section 2.1 we recall 
some results from papers, treating higher current moments or cummulants, es- 
pecially in metallic regime, and the first steps towards 3D criticality. Then we 
add theoretical predictions for any normalized moment and cummulant at 3D 
MIT. In the next Section 2.2 we briefly introduce Nazarov's microscopic model, 
yielding even more information about the moments in metallic regime. In Sec- 
tion 3 we confirm virtually any given analytical formula for current moments 
and cummulants by numerical simulations (except ID and QID). 

2 Model, Numerical Method and Theory 

Let us consider the Anderson model with the Hamiltonian 

H ^^ e„\n >< n| -f ^ \n X n\ + \n' >< n\ (5) 



where the diagonal disorder e„ is taken randomly from < — VK/2,T4^/2 >, i.e. 
it has uniform (box) distribution, and [n, n'] are nearest neighbor sites of a 
d-dimensional cubic lattice. Concerning numerical calculations, the transfer 
matrix method of the calculation of moments and especially the way how to 
deal with larger systems was described in [9] . 

2.1 Limiting values of Fano factor, Cm and current cumu- 
lants 

Lyapunov exponents A„ are related to the eigenvalues 7^ of T as follows: 

r„=.cosh-2(A„) (6) 

It has long been known, that the smallest (most important) A„ exponents are 
distributed equidistantly in large metallic systems, i. e. that their distribution 
function P{\) is constant. In the work [7] this was used to calculate general 
diffusive Cm- Constant density was observed numerically and approved by 
various analytical treatments. Numerical data, however, indicate that in the 
critical regime (3D, 4D) P{\) reaches the origine, and it is not constant. To 
generalize the results, we can propose that at least for smaU A, the P(A) oc A". 
We know already that a = in metallic regime, irrespective of dimension. 
Replacing the trace, i. e. summation of eigenvalues, by an integral over P{\)d\ 
- which works for large enough systems - we get: 

Now let us summarize the metallic regime. The one-third-suppression F = 
1 — C2 = 1/3 was confirmed in [7]. We can also reproduce the Eq. (1), namely 
1 — 3C2 -I- 2C3 = 1/15. The term proposed in [15] is nothing but the fourth 
cumulant: F - 6(C2 - C3) + 6(C3 - C4) = -1/105. And generally for the 
Fano-factor-like quantities we get Cm — Cm+i — Cm/{'^M + 1), with 

^^(^^- (2M-1)! ^^) 

which is yet another practical form of Eq. (22), [7]. An analogy of Eq. (7) was 
proposed already in [9], just with unspecified distribution function. 

In [16] it was shown that a = 1 in 3D at MIT. This was utilized in [13] to 
calculate the limiting value of Fano factor in the critical regime. We can give this 
value exactly: Fc = I — C2(l) = 1/3 4- 1/(6 In 2). In Appendix A the integrals 
of interest are solved analytically. Useful values of Cm{o!) are summarized in 
Table 1. 

The a = 1 is correct not only for -P(A), but also for the distribution of the 
first (smallest) Lyapunov exponent P(Ai), even with the same slope. That is 



why the fits are perfect for all moments - for numerically treated finite systems 
the higher moments depend almost exclusivelly on P(Ai). The possibility of 
a = 2 was proposed in [17] for 4D at MIT; we will come to this point later. 

Let us return to the current cumulants. The generating function Inxaiu) 
from Eq. (3) was calculated for the a = case in [11]: 

Inxo(w) ^Qo dWnl ^ hi = Qo arcsinh^Ve" - 1 (9) 

Jo ^cosh (A) ^ 

where the charge Qo involves constants and physical units, see [11], and it 
is non-universal, depending on model parameters like disorder W. Using Eq. 
(3) they calculated the cumulants <C Q ^ /Qo — 1, ^ Q^ ^ /Qo = 1/3, 
< Q3 > /Qo = 1/15, see Eq.(l), < Q"* ^ /Qo -< Q^ > /Qo = -1/105, see 
above, etc. In the case a = 1 we get a rather formal formula (Appendix A) 

\nxi(u) = Q[J Aln(-£j^ + l)dA=.-2g^y (iln2 - L(t))dt = 

= Q'o In 2 arcsinhVe" - 1 + 2Qo / L{t)dt (10) 

Jo 



where t* = i arcsinhye^^^ and the Lobachevsky function, see [18]: 

L{z) =-£ Hcosix))dx ^^ + ^ + ^ + 0{z^) (11) 

We can calculate the cumulants at 3D criticality, taking again the derivatives 
of Inxi(u) at u ^ 0: 

<Q> = Qo ln2 = Ti (12) 

«Q2» = Q^ (lln2 + i) =Ti-r2 (13) 

«Q'» = Qo (^ln2+^) =ri-3T2 + 2T3 (14) 



« Q* » = Qo (- Y5^ In 2 + — j = Ti - 7T2 + I2T3 - 6r4 (15) 

We added relations in terms of Tm from Eq. (4) . The above-mentioned value 
Fc =< Q^ > / < Q >• 

2.2 Diffusive Tm values from Microscopic theory 

The diffusive contribution to eigenvalues of T was calculated in [14]. Nazarov 
introduced a function dF{(p), and he showed that P(A) is proportional to real 
part of sine/) SF{(f>) at special complex cj) and that 



dF{(P) = Tr 



r 



1 -Tsin^ I 



2(f) 



E' 



(16) 



Here the Cooperon mode energies s^ = ■n'^{n\ +713 + ... + n^), ni = 1, 2, ..., 
712 = 0,1,... , ..., Ud = 0,1,... for simple d-dimensional cube with hard wall 
boundary conditions. The index ni refers to one chosen transfer direction. The 
minor changes for samples with periodic boundary conditions were described e. 
g. in [19]. There is analyzed the 3D case of G = Ti, because the sum on the rhs 
of Eq. (16) diverges for rf > 2 and a cut-off on wave vectors' length is useful. 
Some other details can be found in [20] . 4> is the global phase shift between left 
and right side of the sample and it also plays the role of potential difference [14]. 
We will use it to get the moments of T. Taking derivatives with respect to (fp' 
at zero, details in Appendix B, we get the following linear combinations Rm 



Rr — 



R\ = — T\ ^ 2di 

2 G 

R2 = —72 + -Ti = S — — = 85*2 

i?3 = -T3 + 7^2 - 4^1 = 32^3 (17) 
15 

Ri = -45T4 + 6OT3 - I8T2 + |ti = 5760S'4 

i?5 = -315r5 + 525r4 - 245r3 + —T2 - -Ti = 161280S'5 

3 9 

-14175T6 + 28350r5-17955T4 + 3840r3-186T2 + — Ti = 290304005*6 

-2M 



We introduced the sums Sm = X^s s^'^*^. The way how to calculate these 
sums efficiently in any dimension is also in Appendix B. The first two of these 
equations were derived and the second one compared to numerical simulations 
in [12] and [13]. The values on rhs are summarized in Table 2. 

Recall that for d > 2 the conductivity Ti ^ 00 and for d < 2M the sums 
Sm are finite, see Appendix B. One can divide all the equations by Ti and 
get a homogeneous system for any Cm- Of course the results are the same as 
those from Eq. (22). This means, that Nazarov's formula (16) incorporates the 



Table 1: Selected values of Cm from Eq. (7). Analytical values of Cm(1) can 
be deduced from Appendix A 



a 


C2 


C3 


C4 


C5 


Ce 


C7 


Cs 




1 


2 
3 

0.426217 


8 
15 

0.268839 


16 
35 

0.196084 


128 
315 

0.154259 


256 
693 

0.127120 


1024 
3003 

0.108094 


2048 
6435 

0.0940195 



diffusive case a = of Eq. (7), though the validity of the last goes beyond the 
mentioned restrictions on dimension (or M). Therefore it is usefuU to calculate 
1 — C2 by numerical simulations at 4D and show that the 1/3 suppression holds, 
see Fig. 3 left. If we would like to derive this analytically also for d > 2M, 
the cut-off K could be introduced together with finite sums Sm{K), where the 
summation is performed only over finite set of {rii, ...jfi^;} such that the wave 
vector length s < K. Then we expect lim^^oo S m{K) / Si{K) — for all M > 2 
and d > 2. 

To proceed with comparing Eq. (16) and the attitude from previous subsec- 
tion, we would need the - dependence of the Lyapunov exponents distribution 
P(A, (j)). We can conjecture that it is independent of A, i. c. P(A, 0) = G(0), and 
G{0) = Ti. Inserting this into the first expression for the generating function in 
Eq. (16) and replacing the trace as in Eq. (7), we get: 



SF{<P) = G(0) 



cosh(A)-2(iA 
1 — cosh(A)^2 31^2 



G(</>) dX 



cosh(A) 



sm0 

(18) 

which is consistent with Eq. (17) at (f> ^ 0, and comparison to the rhs. of Eq. 
(16) yields Gicj)). Nazarov [14] introduced yet another quantity I{(f>) = (p G and 
he gave a formula for the Fano factor in terms of derivatives of /. We are now 
able to evaluate them: 



1 



2 r"(o) 
~p(or 



1 



12^2 

Si 



(19) 



We can see again, that the one-third suppression is given by the infinite de- 
nominator 5*1 rather than by zero numerator (constant, (j) independent G), as 
previously assumed. 

It also is interesting, that the ^2 value appears in the formula for UCF, e. 
g- [20], [21]: 



{G^)-{Gf^l2S2^-{S 



3 / 



(G) - 2^2 



(20) 



Table 2: The sums Sm with appropriate factors in dimensions 2, 3 and 4 
d 8 S'2 32 5*3 5760 ^4 161280 5*5 29030400 Sq 

2 0.123742 0.0387682 0.649759 1.77893 31.9117 

3 0.209369 0.0458759 0.699664 1.84186 32.4521 

4 00 0.0573425 0.763394 1.91503 33.0469 
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Figure 1: Fano factor 1 — C2 = F and C3 as a function of system size. The 
description of disorders applies to both figures. Dashed Hues are theoretical 
diffusive values from Table 1. 

3 Numerical simulations 

We calculated the mean values of Tm and Cm, M ~ 1,2,..., 8 in dimensions 
rf = 2,3 and 4 for various disorders and system lengths L. Typical ensembles 
were of 10^ elements for smaller systems and 3 x 10^ for the greatest ones. 



3.1 Fano factor and Cm 

d = 2 is lower critical dimension for Anderson model and the case with or- 
thogonal symmetry has still no phase transition. Any positive disorder forces 
that conductivity tends to zero for very large systems. The obvious inequalities 
< Tm+1 < Tm < ■■■ <Ti — G do not mean, that Cm should necessarily tend 
to zero. Nevertheless this can be expected for extremely large samples, as they 
approach the localized (Poissonian) regime. From Fig. 1 we see, that C2, C3 are 
monotonic - decaying - functions of L and that the theoretical metallic values 
are close to the plateau-like region, where the L-dependence is weak. But there 
is neither any special disorder, nor an asymptotic value, that would be reached. 
The C4, ...,Cs,{L) dependences look very similarly, just as if they were scaled 
down to lower values. The Cm{^/L) dependences for H^ = 3 and M = 2, ...,8 
in 2D were pubHshed already in [9] (and compared to Eq. (22) in [7]), together 
with 3D data for W = 10 and Wc = 16.5. 

The situation changes substantially in higher dimensions d > 2. The metal- 
lic limit becomes attractive for disorders < M^ < Wc- The MIT is manifested 
by a constant line at W — Wc, which repels neighboring lines on both sides, 
i. e. also for W > Wc, slowly approaching Poissonian limit. This is why we 
think that any Cm is a good order parameter, the instability of the critical 
point follows from the scaling hypothesis. In 3D both the metallic and critical 
values agree very well with theoretical predictions of Eq. (7) for a = and 1, 
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Figure 2: C3 and C4 as a function of system size in 3D. The description of 
disorders applies to both figures. Dashed fines are theoretical values from Table 
1. 
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Figure 3: Fano factor F = 1 — C2 and C3 as a function of system size in 4D. 
Dashed lines are theoretical metallic values from Table 1. The description of 
disorders applies to both figures. 



respectively, see Fig. 2. Fano factor F{L) has already been published in [13]. 
Concerning 4D, metallic limits are still the same, as expected, see Fig. 3. But 
the MIT values are not so easy to understand. We tried several conjectured 
P(A), e. g. with non-integer a, or a polynomial. None of them could satisfy 
all C'm values with low number of fitting parameters, though polynomials of 
the type with Pc(A) = biX + 62A^ + b^X^ can fit the moments within several 
percent. Most probably we have to do with some nontrivial function, too com- 
plicated to be found numerically. Not to speak about the possibility of different 
-Pc(Ai). Of course, in 4D we can still have huge finite size effects, spoiling the 
precision of Cm values at MIT. We leave the 4D Pc{X) question open for future 
investigations. 
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Figure 4: Various Rm quantities as a function of conductivity in 2D. Open 
symbols are R2 and R4, full ones R^ and R^. Dashed lines are theoretical 
values from Table 2. 



3.2 Diffusive Rm combinations of Tm 

We calculated also the averaged Rm quantities in order to check Eq. (9) numer- 
ically. It is usefuU to plot them as a function of conductivity G, rather than of 
L. Fig. 4 left shows, that in 2D the R2 only roughly approaches the predicted 
value, at larger G (smaller W) the ballistic regime disturbs the convergence. 
i?2(G) was shown already in [12]. In order to judge the influence of ballistic, we 
plot the data also for lower W and L in Fig. 5. Without any disorder {W — 0) 
all 7^ = 1 and thus all Tm = Ti. So we get the ballistic asymptotics, e. g. 
— G/3 and — 2G/15 for R2 and i?3, resp. One can see that the data for low W 
and L are close to the ballistic asymptotics and with rising L they approach the 
metallic limit. Surprisingly, for the higher indexed i?3, _R4, i?5, ... this problem 
becomes less sensitive. Especially if we limit ourselves only to data for the sev- 
eral largest L, their behaviour is almost perfect. Fig. 4, within few percent from 
theoretical values. Table 2. To overcome the ballistic influence for R2 fully, we 
would need even greater L than available. 

In 3D there are no signs of disturbing ballistics even for R2 and our data 
are in better accordance with predicted value than those in [13]. Again the Eq. 
(16) is valid without any doubt, see Fig. 6. 

The i?2 in 4D is expected to diverge. This is not easy to show numerically, 
as available system sizes are rather small. Fig. 7 left, and ballistic influence 
is large again, similarly to 2D. Retrospectively we might even say, that in 2D 
the i?2 alone was insufficient to confirm any statement. Back in 4D we can 
just state, that simulated data are in no contradiction with possible logarithmic 
divergence R2 oc ln(G) ex ln(L), which we would expect from analogy with 4D 
UCF, see [20]. On contrary, i?3,i?4... converge again unexpectedly well to the 
predicted limits. Fig. 7. 

Summarizing, we regard Nazarov's formula (9) as confirmed by numerical 
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Figure 5: R2 and R^ in 2D including smaller W and L. Horizontal dashed lines 
are metallic theoretical values from Tabic 2. Long - dashed lines are ballistic 
limits. The description of disorders applies to both figures. 



simulations in Anderson model at least for dimensions up to 4. 

4 Conclusions 

We have shown all normalized higher current moments Cm can be calculated 
exactly in the large system size limit at least in next regimes: the diffusive one in 
virtually any dimension [7], and the critical one (MIT) in 3D. We have confirmed 
the predicted values by numerical simulations. The MIT in 4D remains an open 
question. Further we saw, that Nazarov's formula (16) yields the same diffusive 
Cm values as our Eq. (22) and also special linear combinations of the Tm 
moments converging to nontrivial constants. We have convincingly confirmed 
the validity of Eq. (16) for Anderson model in metallic regime by numerical 
simulations. 



5 Appendix A 

The aim is to find analytical solution of the integral 

^M(a) = 



cosh (x) 



dx 



(21) 



appearing in the Eq. (5) for Cm — Im /Ii- Special cases known from litera- 
ture [18] are a — Q (yielding directly the Cm(0) as /i(0) = 1, see [7]): 



/m(0) = Cm(0) = 4^'^-^B(M,M) - -B(1/2,M) 



r(i/2)r(M) 

2r(M + 1/2) 



(22) 
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Figure 6: Various Rm quantities as a function of conductivity in 3D. Open 
symbols are R2 and i?4, full ones R^ and Re,. Dashed lines are theoretical 
values from Table 2. The description of disorders applies to both figures. 



and M=l: 



h{a) = 2^-°'{l - 2i"")r(a + l)C(a) if a 7^ 1 



(23) 



We introduced the special functions F, B and C, i. e. the Gamma, Beta and 
Rieman's Zeta function, respectively, [18]. Further we use the polylogarithmic 
function LiQ,(x): 



Li«(a;)-I]f^ 



k=l 



more specifically its value at x — —1: 



(24) 



Li„(-1) = -(l-2i-")C(a) if a 7^1 
- In 2 if a = 1 

Without going into details, we found the following exact formulas: 



(25) 
(26) 



h{a) 
Ii{a) 



-2i-"r(a + l)Li„(-l) (27) 

^ "F(a + l)[Li„_2(-l)-Li„(-l)] (28) 

T{a + l)[Li„_4(-l) - 5Li„_2(-l) + 4Li„(-l)] (29) 



3! 

25-Q 



5! 

27-a 



V{a + l)[Li„_6(-l) - 14Li„_4(-l) + 49Li„_2(-l) 



7! 
-36Li„(-l 



(30) 
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Figure 7: Various Rm quantities as a function of conductivity in 4D. Open 
symbols are R2 and R4, full ones R3 and R^. Dashed lines are theoretical 
values from Table 2. The description of disorders applies to both figures. 

/5(a) = -^r(a + l)[Li„_8(-l)-30Li„_6(-l) + 273Li„_4(-l) 

-820Li„_2(-l) + 576Li„(-l)] (31) 

nil — Q: 

/6(a) = ^Y^r(a + l)[Li„_io(-l)-55Li„_8(-l) + 1023Li„_6(-l) 

-7645Li„_4(-l) + 21076Li„_2(-l) - 14400Li„(-l)] (32) 

It is worth mention that LiQ,(— 1) and C(q^) have rational or zero value for a 
negative integer or zero. In particular, /i(l) = In 2, hi^) = — ^ + | ln2, 13(1) = 
— ii + :^ ln2, 14(1) = —j^ + M ^'^^' ^*'^- These rationals (the term at In 2 is 
/m(0)) can be found much easier with the help of the generator [18]: 

, , , r°° X dx t\n2- L(t) 

lait) = / ; ;- = 4 . ^ (33) 

Jo cosh X + cos 2t sm 2t 

One just takes consecutive odd derivatives aX t — 7r/4. It is interesting to 
compare this generator to Eq. (18). Integrating both sides of Eq. (33) w. r. to 
y = — sin^ t yields Eq. (10). 

6 Appendix B 

In order to get relations for Tm, we will perform n-fold derivation (n = 0, 1, ..., M- 
1) of 6F{(f>) with respect to 0^ at = 0, on both terms from Eq. (16): 

hm SF = Ti = -2Si (34) 

lim^5F=limi-A5F = Il = -2S2-^ (35) 

0^0 d(j)^ <t>^0 2(j)d4: 4 3 ^ ' 
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etc. We solve this linear system of equations so that we remove almost all the 
sums Sn, except of the one with the highest index. Thus we get the set of Eq. 
(17). We choose this one from the plenty of possibilities, because Sm diverges 
for d > 2M, like S2 and R2 at 4D, Fig. 7. But S3 is finite at 4D (even 5D) and 
it can be compared to numerical results. The case M = 2 was already treated 
this way in [12] and [13] and the sums S2{d) were calculated for d = 2,3. 

Now let us show how to calculate the d-fold sums Smid) easily with much 
better precision than by direct summation. The statement is: 

^M(d) = ^^^g^y^ j^^—^[M0,e--n~W3{0,e--y) + lf-^dy (37) 
Here i^s is Jacobi elliptic function, simplified to: 

M0,,)= ±^'' = j=± -p(-^) (38) 

with q — cxp(— Try). The proof of Eq. (37) consists in using the first expansion 
of elliptic functions and performing per partes (M — l)-times. This time the case 
M = 2 was proposed in [20] , dealing with universal conductance fluctuations in 
arbitrary dimension. 

Finally let us derive the convergence condition for SM{d). The values of 
the expression under integral (37) at small y are crucial. We use the second 
expansion -d^lO^q) « l/^/y (case fc = 0) and the leading term is proportional 
to yAf-i-d/2_ r^Yic turning point of convergence is y~^, i. e. M = d/2, with 
logarithmic type of divergence [20] . 
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